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Abstract 


A significant excess of the stellar mass density at high redshift has been discovered from the early data release of 
James Webb Space Telescope (JWST), and it may require a high star formation efficiency. However, this will lead 
to large number density of ionizing photons in the epoch of reionization (EoR), so that the reionization history will 
be changed, which can arise tension with the current EoR observations. Warm dark matter (WDM), via the free 
streaming effect, can suppress the formation of small-scale structure as well as low-mass galaxies. This provides an 
effective way to decrease the ionizing photons when considering a large star formation efficiency in high-z massive 
galaxies without altering the cosmic reionization history. On the other hand, the constraints on the properties of 
WDM can be derived from the JWST observations. In this work, we study WDM as a possible solution to 
reconcile the JWST stellar mass density of high-z massive galaxies and reionization history. We find that, the 
JWST high-z comoving cumulative stellar mass density alone has no significant preference for either CDM or 
WDM model. But using the observational data of other stellar mass density measurements and reionization history, 
we obtain that the WDM particle mass with mw = 0.5170: 12 keV and star formation efficiency parameter 


ris > 0.39 in 2c confidence level can match both the JWST high-z comoving cumulative stellar mass density and 


the reionization history. 
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1. Introduction 


High redshift (high-z) galaxies play an important role in the 
Epoch of Reionization (EoR), which is considered to be the last 
phase transition of the Intergalactic Medium (IGM) in the 
cosmic history. Studying the high-z galaxies allows us to have 
a deeper understanding on the formation and evolution of early 
galaxies and the cosmic reionization history. The launch of the 
James Webb Space Telescope (JWST) opens a new era to study 
the high redshift object in the Universe (Finkelstein et al. 
2022, 2023). The JWST early observational data release 
provides valuable information about the properties of high-z 
galaxies, which may challenge the current galaxy formation 
theory under the ACDM model (Menci et al. 2022; Naidu et al. 
2022a, 2022b; Boylan-Kolchin 2023; Mason et al. 2023; 
Mirocha & Furlanetto 2023). 

In particular, Labbé et al. (2023) may find a serious tension 
with the standard galaxy formation theory. Using the JWST 
Cosmic Evolution Early Release Science sample (detected by 
the JWST/NIRCam instrument), they identified six galaxies 
with stellar masses My = 10? ~ 10! 1 ! M. at 7.4 «& z « 9.1. 
The corresponding comoving cumulative stellar mass density is 


about 20 times higher at z~8 and about three orders of 
magnitude higher at z~ 9 than the prediction from the star 
formation theory in standard ACDM cosmology, which is 
based on previous observations (Stark et al. 2013; Song et al. 
2016; Bhatawdekar et al. 2019; Kikuchihara et al. 2020; 
Stefanon et al. 2021). This huge abundance excess may due to 
issues of galaxy selection, measurements of galaxy stellar mass 
and redshift, dust extinction, and sample variance (Endsley 
et al. 2022; Ferrara et al. 2022; Adams et al. 2023; Ziparo et al. 
2023). But if this result is confirmed by future spectroscopic 
observation, for example, the follow-up observation by the 
JWST/NIRSpec, it will be a huge challenge to the ACDM 
model. There are many efforts trying to solve this tension, for 
example, the primordial non-Gaussianity in the initial condi- 
tions of cosmological perturbations (Biagetti et al. 2023); the 
rapidly accelerating primordial black holes (Liu & 
Bromm 2022; Yuan et al. 2023); the Fuzzy Dark Matter 
(FDM) (Gong et al. 2023); axion miniclusters or primordial 
black holes (Hiitsi et al. 2023); the cosmic string (Jiao et al. 
2023); a gradual transition in the stellar Initial Mass Function 
(Trinca et al. 2023); modified ACDM power spectrum 
(Padmanabhan & Loeb 2023; Parashari & Laha 2023). 
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A direct way to explain the observation of Labbé et al. 
(2023) in theory is enhancing the star formation efficiency f... 
Some previous observations prefer a low star formation 
efficiency (fx « 0.1). But since f, depends on complicated 
astrophysical processes which remain unclear, it is theoretically 
possible to have a larger f;. at high redshifts. On the other hand, 
although assuming large fẹą at high redshifts may help to solve 
or reduce the tension between the current model and the JWST 
Observation, it may raise new problems. If we have higher star- 
forming efficiency at high-z, the total number of ionizing 
photons in the EoR will significantly increase, and then the 
EoR may end much earlier which could be in tension with the 
EoR history observations (Hinshaw et al. 2013; Planck 
Collaboration et al. 2020). A possible way to resolve this 
problem is to enhance the formation rate of massive galaxies, 
but suppress the formation of low mass galaxies at the same 
time. A mass-dependent star formation efficiency f,.(M) can 
suppress the star formation in small halos, but may still not be 
enough to reduce the huge tension. 

In Gong et al. (2023), they propose to use FDM to solve this 
problem. The FDM is expected to have extremely low mass of 
10 7" < m, X 10? eV, and it can suppress the abundance of 
low-mass haloes via the effective quantum pressure that arises 
from its galaxy-size de Broglie wavelength, which can 
successfully explain both JWST and EoR observations. In this 
work, we study warm dark matter (WDM) as an alternative 
solution, and explore the implications on the WDM properties 
with the JWST data. In the WDM scenario, WDM particles 
with mass of a few keV are much lighter than the standard 
CDM particles, allowing them to remain relativistic for longer 
in the early universe and to retain a non-negligible thermal 
velocity dispersion. This thermal velocity dispersion will 
cause a so-called free streaming effect, making WDM 
particles escape out of the high density region and resulting 
in a suppression of the structure growth on small scales 
(Blumenthal et al. 1982; Bode et al. 2001; de Vega & 
Sanchez 2012). Then the abundance of low-mass haloes as well 
as galaxies will be significantly suppressed, and this highlights 
us to consider WDM as a solution to explain the JWST 
observation (Labbé et al. 2023). In this work, we assume 
a flat universe with Qm = 0.3153, Q, = 0.0493, h = 0.6736, 
0g = 0.811, n, = 0.9649 (Planck Collaboration et al. 2020). 


2. Model 
2.1. WDM Halo Mass Function 


In order to analyze the impact on the halo abundance by 
WDM, we need to calculate the halo mass function. A 
conventional choice for the halo mass function is obtained from 
the ellipsoidal collapse model (Press & Schechter 1974; Sheth 
& Tormen 1999; Sheth et al. 2001; Cooray & Sheth 2002), 
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which can be expressed as 
n(M)dM — af dy, (1) 


where M is the halo mass, p is the mean comoving matter 
density, and 


fv) — AD + (qare, (2) 


Here p = 0.3, g = 0.707, A = 0.2161 are the parameters 
derived from simulations, and (M) = 6./o(M). 

In the CDM case, we assume the critical overdensity is 
independent of the mass scale, and 6, is the critical overdensity 
barrier for collapse, which can be derived analytically and 
given by 6. e 1.686 for the spherical collapse of CDM. On the 
other hand, in the WDM case, we expect that the free streaming 
effect will allow the WDM particles escape out of the 
collapsing region, and as a result, the collapse of WDM will 
become more difficult on small scales. Barkana et al. (2001) 
studied the collapse thresholds for WDM by performing 
spherical collapse simulations. They find that, when the halo 
mass is below a characteristic mass scale, the threshold for 
collapse will increase rapidly with decreasing mass. Benson 
et al. (2013) provided a fitting function to describe the mass 
dependence of 6, as 


6° 
6-(M, z) = GM) —, 3 
(M, z) = GM) DO (3) 
where 
0.31687 
G(M) = (1 — hG))exp | 03095 | 
0.04 
T0 exp(2.3x) d 


Here x = log(M/M)), and M; is the effective Jeans mass of the 
WDM as defined by Barkana et al. (2001): 


. 2 
Mj = 3.06 x DE a y (Se 


3000 0.15 
=i —4 
x (&) mw ) Mo, (5) 
1.5 ] keV 
where Zeg = 3600( 2a — 1 is the redshift of matter-radiation 


equality, gw = 1.5 is the effective number of degrees of freedom 
for the spin -5 fermion, and h(x) = 1/[1 + exp(10x + 24)] is 
an auxiliary fitting function. c(M) can be calculated as the 
variance of the linear matter overdensity field when smoothed on 
a comoving scale R: 


o6 kW. 
o? (R) = f an() Plin(k) W2(kR)d Ink, (6) 
0 2m 


where W (x) — (3/x?)(sinx — x cosx) is the Fourier transform 
of a spherical top-hat filter window function. We can rewrite 
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o(R) in terms of mass by relate the comoving scale R with the 
the mass enclosed within this scale as M = (4/3) 7R?p. P'"(k) 
is the linear matter power spectrum, and it can represent the 
CDM and WDM linear power spectrum in this study. 

Since WDM can suppress the matter clustering on small 
scales via the free streaming effect, the WDM linear matter 
power can be calculated by Bode et al. (2001): 


Pipa (k, z) = Paba(k, z) Tw (k). (7) 


The CDM linear power spectrum Pcpm(k, z) is obtained by 
CAMB (Lewis et al. 2000), and the transfer function Tw(K) is 
assumed to be redshift independent, and it is given by 


Tw(k) = [1 + (oy? ] 7n. (8) 


Here the fitting parameters are 
0.074(mw/1 keV)--(0.7/h) Mpc. 

In Figure 1 we show the halo mass functions for both WDM 
and CDM cases at redshift z= 8, and we can find that the 
abundance of the low-mass halos are strongly suppressed by 
the free streaming of WDM. So the contribution of ionizing 
photons from low-mass galaxies will also be suppressed, and 
we are allowed to have a larger number density for massive 
galaxies without altering the reionization history. 


p-2112, a= 


2.2. Stellar Mass Density and Reionization History 


After obtaining the halo mass function, the comoving 
cumulative halo mass density with halo mass greater than M 
can be estimated by: 


p(>M, z) = J dM'M'n(M', z). (9) 


By considering the M,—M relation My = (%»/9m M = 
c M, the cumulative stellar mass density with stellar mass larger 
than M. can be written as : 


DX >Mx, z) = ep M«fe, z). (10) 


The star formation efficiency f, can be considered as a mass- 
dependent quantity. In this work, we assume a double power- 
law form of f,(M), and following Mirocha et al. (2017) we 
have 


2n 


M Mo M Qhi ? 


where R is the star formation efficiency at its peak mass Mp, 
and oj, and api describe the power-law index at low and high 
masses. In our work, we adopt M, = 2.8 x 10!! Mo, Ojo = 0.49 
and ay; = —0.61. The parameter i, will be adjusted to match 
the JWST data in this work. 

The star formation efficiency is tightly related to the 
emission of ionizing photons, and thus can change the cosmic 
reionization history. Usually, we can explore the reionization 
history by investigating the redshift evolution of the hydrogen 


fM) = (11) 
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Figure 1. The halo mass functions at redshift z — 8 for different WDM masses 
(solid lines). The CDM case is also shown for comparison (black dashed line). 
We can see that WDM can significantly suppress the formation of low-mass 
halos, and the suppression mass scale increases with the WDM particle mass 
decreases. 


volume filing factor Qyy(z) by solving the differential 
equation as below (Madau et al. 1999; Wyithe & Loeb 2003): 


d lion = 
Mun — 6, Mem — Gyu(aw (Tu Aal + Dre (12) 
dt ny 


Here we set the escape fraction fese = 0.1 (Sun et al. 2021), fig 
is the mean number density of hydrogen (both neutral and 
ionized) atoms today, Cg n = 3.0 is the clumping factor of the 
ionized gas (Kaurov & Gnedin 2015) which is assumed to be 
redshift-independent here, og is the Case B recombination 
coefficient, and Tyn is the kinetic temperature. We set 
TH n= 2 x 10* K as a constant (Robertson et al. 2015), and 
then the recombination coefficient ap = 2.5 x 107} em? & *. 
The total ionization fraction x, then can be estimated as below, 
by assuming that the helium has the same first stage ionization 
fraction as hydrogen, and we have 


ve ox! m e), (13) 


where Ype = 0.25 is the helium element abundance. The 
ionizing photons emission rate per unit comoving volume ion 
can be estimated by 


Oy 1 oo 
E" > Z) 14 
Ren Ja. dMf,(M) Mn(M, z) (14) 


where Nion 4000 is the total ionizing photon number that a 
stellar baryon can product throughout its lifetime for typical 


Nion = N; 
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The comoving cumulative stellar 
mass density at z = 8 
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The comoving cumulative stellar 
mass density at z = 9 
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Figure 2. The comoving cumulative stellar mass density for different WDM masses and the CDM case at z = 8 (left panel) and z = 9 (right panel) with fe = 0.5. We 
also show the range of the results of my = 0.6 keV with fe = 0.1-1.0 with the orange hatched region. The JWST results from Labbé et al. (2023) at z ~ 8 and 9 are 
shown by the gray shaded regions. We also compare our result with some previous measurements (Stark et al. 2013; Song et al. 2016; Bhatawdekar et al. 2019; 
Kikuchihara et al. 2020; Stefanon et al. 2021) at low stellar mass region shown as data points with error bars. 


Population II galaxies (Leitherer et al. 1999; Vazquez & 
Leitherer 2005; Leitherer et al. 2010, 2014). tsp is the star 
formation timescale, which equals 10% of the Hubble time at 
redshift z (Wyithe & Loeb 2006; Lidz et al. 2011). Mnin stands 
for the minimum halo mass corresponding to a virial temperature 
Tyi» and halos above this mass can sustain effective cooling via 
the Lya transition (Barkana & Loeb 2001): 


l 2/3 
Ti = 1.98 x (£) — Monin 
0.6 J^ 1051 Ms 


1/3 
: | Qm Ac | (L5). 
Qm(z) 187? 10 

where u = 0.61 is the mean molecular weight, Qm(z) is energy 

density fraction of matter component at redshift z, and 

A,—18m^--82d —39d^ with d=Qy(z)—1. We take 

Tirs 10* K to obtain the corresponding Mmin at each redshift. 

In Equation (14), we can find that the shape of £.(M)n(M, z) 

mainly determines the contribution of ionizing photons from 

galaxies with different masses. In the CDM model, the low- 

mass galaxies will make a great contribution, but in WDM 

case, since the formation of low-mass galaxies is suppressed by 

WDM free streaming effect, the ionizing photons will be 
dominated by massive galaxies. 

In addition, the optical depth of the cosmic microwave 


background (CMB) scattering is also a good probe for 
characterizing the reionization history, and it can be estimated by 


cdz 
“(+ 2H(2)’ 


(15) 


qum L orñyg(l + z)x (16) 


where or= 6.65 x 10 7 cm? 


cross-section. 


is the Thompson scattering 


3. Result 


In Figure 2, we present the comoving cumulative stellar mass 
density at redshift z= 8 and z = 9 for different WDM masses and 
the CDM case. We find that both WDM and CDM could 
consistent with JWST data when the star formation efficiency 
parameter JE ~ 0.5. But if we compare the result with the 
previous measurements (e.g., Hubble Ultra-Deep Field, Spitzer/ 
Infrared Array Camera, Keck Observatory, and so on) at low 
stellar mass region (Stark et al. 2013; Song et al. 2016; 
Bhatawdekar et al. 2019; Kikuchihara et al. 2020; Stefanon 
et al. 2021), we can see that the CDM model will significantly 
over predict the stellar mass density at this region, since too many 
low-mass galaxies will form under such high star formation 
efficiency. On the other hand, the advantage of the WDM model 
is that it will strongly suppress the formation of low-mass 
galaxies, and is able to match the measurements of both low-mass 
and high-mass regions. As a reference, we also show the result by 
varying the star formation efficiency parameter Li — 0.1-1.0 
with my — 0.6 keV shown as a hatched orange region. 

As we discuss in Section 2.2, we can also check our result by 
investigating the impact of our model to the cosmic reionization 
history. We calculate the optical depth 7 as a function of redshift 
for the WDM and CDM models, as shown in the left panel of 
Figure 3, and compare it with the 9 yr WAMP (Hinshaw et al. 
2013) and Planck 2018 (Planck Collaboration et al. 2020) results. 
We can see that the model with my=0.6keV is in good 
agreement with Planck 2018 measurements. 

In the right panel of Figure 3, we present our calculation result 
of the neutral hydrogen fraction characterized by 1 — Qu n as a 
function of redshift. The data of observational constraints on the 
neutral hydrogen fraction from different methods are also 
shown for comparison, such as Lya forest dark fraction (McGreer 
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Optical Depth with f} =0.5 
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Figure 3. Left panel: the optical depth 7 as a function of redshift for different WDM masses and the CDM cases with R = 0.5. The results of mw = 0.6 keV with 
fe = 0.1-1.0 is shown as the orange hatched region. We compare our results with the 1c (68% C.L.) measurements from WMAP (Hinshaw et al. 2013) and Planck 
(Planck Collaboration et al. 2020), which are shown as gray shaded bands. Right panel: the neutral hydrogen fraction denoted by 1 — Qy n as a function of redshift. 
The data points are the results from previous observations (Totani et al. 2006; Ota et al. 2008; Dijkstra et al. 2011; Stark et al. 2011; McGreer et al. 2015; Sobacchi & 
Mesinger 2015; Davies et al. 2018; Mason et al. 2018, 2019; Greig et al. 2019; Hoag et al. 2019). 


et al. 2015), Lya emitters (Ota et al. 2008), Lya equivalent width 
(EW) distribution (Mason et al. 2018, 2019; Hoag et al. 2019), Lya 
galaxy clustering (Sobacchi & Mesinger 2015), Lya emitter 
fraction (Dijkstra et al. 2011; Stark et al. 2011), QSO damping 
wing (Davies et al. 2018; Greig et al. 2019), and GRB damping 
wing absorption (Totani et al. 2006). We can see that, the WDM 
model with mw =0.6keV and fe = 0.1-1.0 (shown by the 
orange hatched region) is in good agreement with those 
observational data. So our result indicates that the model with 
mw ~0.6keV can match both the high redshift comoving 
cumulative stellar mass density and the reionization history. 

We also sample the posterior distribution of the model 
parameters using the Markov Chain Monte Carlo (MCMC) 
method with the Cobaya package (Torrado & 
Lewis 2019, 2021). The flat priors are used for the two model 
parameters with 0.1 < mw/keV < 3.0 and 0 < Ps < 1. Since 
there is no effective constraint with JWST data only, we 
constrain the parameters using the cumulative stellar mass 
density data from both JWST and previous measurements at 
z=8 and 9, and also consider the reionization history data, 
including the optical depth measured by Planck and the neutral 
hydrogen fraction measured by various observations (see the 
data in Figures 2 and 3). 

We present the constraint result in Figure 4. As shown by the 
red contours, when only considering the cumulative stellar 
mass density data from JWST and previous measurements, it 
favors a large my, and a small E which is similar as the CDM 


case, or a small mẹ and a large ro The 1D probability 
distribution functions of the two parameters are shown in red 


curves, and we have m, = 0.511449 (1o) and 0.51277 


(2c) keV and f = 0.147399 (1c) and 0.14173 (20). When 
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Figure 4. The MCMC constraint result on my and fe . Red contours (68.396 
and 95.5% C.L.) and curves show the constraint from the high-z comoving 
cumulative stellar mass density data given by JWST and previous measure- 
ments at z — 8 and 9 shown in Figure 2, and the blue contours and curves show 
the joint constraint result using both the high-z galaxy stellar mass density data 
and reionization history data (as shown in the left and right panels of Figure 3). 


including the data from reionization history, we find that the 
constraint on my, is significantly improved with my = 0.514813 
(10) and 0.517072 (2c) keV, and only large R > 0.51 (1o) and 


70.39 (20) is favored in this case. 
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We notice that, this result may have some tension compared to 
previous studies. For example, using the UV luminosity functions 
of high redshift galaxies, a lower limit of my = 2.9 keV can be 
obtained (Schultz et al. 2014; Menci et al. 2016; Corasaniti et al. 
2017; Rudakovskyi et al. 2021). A similar result can be derived 
using the global 21 cm signal detected by Experiment to Detect 
the Global Epoch of Reionization Signature observation 
(Bowman et al. 2018), where the WDM models with mw < 
3 keV can be ruled out (Chatterjee et al. 2019). We find that, in 
these works, a low star formation efficiency (e.g., fe < 0.1) is 
always assumed, which is derived from the observations of 
Hubble Space Telescope (HST) and other telescopes, and this 
may be different from the result implied by the JWST data given 
by Labbé et al. (2023), which may favor a higher star formation 
efficiency with f > 0.1. In Figures 1-3, we plot the results of 
the WDM model with my — 3 keV for comparison. We can see 
that the WDM model with my = 3 keV has a suppression mass 
scale ~10’M., and obviously cannot fit the data of reionization 
history with a large star formation efficiency (f. » — 0.5). A small 


star formation efficiency with A « 0.1 is needed in this case to 
match the reionization data, but then the JWST data shown in 
Figure 2 cannot be matched, just like the CDM case. Besides, 
Enzi et al. (2021) derived a lower limit of my > 6.048 keV by 
combining the strong gravitational lensing, the Lyo forest and the 
abundance of the Milky Way satellites, which depends on the 
low-mass dark halo mass function and is obviously a challenge to 
our model. In our low WDM mass case, the formation of the dark 
halos below 10? ~ 10° Mo are strongly suppressed, and it may 
be challenging to explain the existence of low-mass galaxies in 
the Milky Way (McConnachie 2012; Newton et al. 2018). Hence, 
it seems that the inconsistency of the WDM particle mass 
between our result and the previous studies probably comes from 
the inconsistency between the data of JWST and other previous 
observations. 

We should note that there are indeed some arguments on the 
data given by Labbé et al. (2023). For instance, Prada et al. 
(2023) found that their galaxy formation model can explain the 
UV luminosity function measured by JWST (Naidu et al. 
2022b; Donnan et al. 2023; Harikane et al. 2023) and HST 
(Oesch et al. 2018; Bouwens et al. 2021; Kauffmann et al. 
2022), and their prediction on star formation rate and stellar 
mass at z~ 8.5 has good match with the observation of the 
spectroscopically confirmed galaxies (Bouwens et al. 2022; 
Fujimoto et al. 2023; Heintz et al. 2023; Williams et al. 2023). 
So they claim that the stellar mass-to-light ratio during early 
epochs could not have reached such high values reported by 
Labbé et al. (2023). Chen et al. (2023) discussed three sources 
of uncertainties in counting massive galaxies at high-z, 
including cosmic variance, error in stellar mass estimate, 
and contribution by backsplash. They found that each of 
them can significantly affect the estimation of stellar mass 
density at high-z. Besides, photometric redshift estimation 
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(Adams et al. 2023) and dust extinction (Ferrara et al. 2022; 
Ziparo et al. 2023) are also important factors that need to be 
carefully calibrated in the data analysis. So follow-up 
measurements by spectroscopic observations (especially by 
the JWST/NIRSpec) are probably necessary for further 
confirmation of these data in the future. 


4. Conclusion 


Unexpected high stellar mass densities of massive galaxies at 
z — 8 and 9 have been found by the early JWST observations, 
which may indicate a high star formation efficiency and can 
have tension with the current measurements of the cosmic 
reionization history. We propose that the WDM model can be a 
solution and explore the implications on the properties of 
WDM particles. Since WDM can suppress the formation of 
low-mass halos by free streaming effect, the number of small 
galaxies and hence the ionizing photons will be effectively 
decreased. Therefore, in this scenario, the reionization history 
would not be changed with a high star formation efficiency for 
high-z massive galaxies. 

After comparing the predictions of our WDM model to the 
JWST (at high stellar masses) and other high-z galaxy 
observational data (at low stellar masses), the CMB optical 
depth measurement from Planck, and the neutral hydrogen 
fraction of IGM measurement from a bunch of different 
methods, we employ the MCMC method to fit those 
observational data, and find that the WDM model with 
mw = 0.51032 keV with x > 0.39 is in good agreement 
with all these observational data in 95.596 C.L. 

We should note that this result still has large uncertainties 
and may change with the adopted models and observational 
data. On one hand, there still may be large uncertainties in our 
current theoretical models and parameters about the WDM, 
optical depth, and reionization history, which can significantly 
change the derived WDM particle mass. On the other hand, the 
early release data of JWST need to be further confirmed by 
follow up spectroscopic observations, especially with JWST/ 
NIRCam. If we can obtain higher quality data from future 
observations of JWST, we are supposed to put a more reliable 
and tighter constraint on WDM particle mass, and have a better 
understanding on the nature of dark matter and the galaxy 
formation process. 
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